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The article addresses morphological approaches to design of modular systems. The following methods 
are briefly described: (i) basic version of morphological analysis (MA), (ii) modification of MA as method 
of closeness to ideal point(s), (iii) reducing of MA to linear programming, (iv) multiple choice problem, (v) 
quadratic assignment problem, (vi) Pareto-based MA (i.e., revelation of Pareto-efficient solutions), (vii) 
Hierarchical Morphological Multicriteria Design (HMMD) approach, and (viii) Hierarchical Morphological 
Multicriteria Design (HMMD) approach based on fuzzy estimates. The above-mentioned methods are 
illustrated by schemes, models, and illustrative examples. An additional realistic example (design of GSM 
network) is presented to illustrate main considered methods. 

Keywords: System design, Morphological design, Modular systems, Configuration, Composition, Syn- 
thesis, Combinatorial optimization, Decision making 



1. Introduction 

Morphological analysis (MA) was firstly suggested by F. Zwicky in 1943 for design of aerospace systems. 
Morphological analysis is a well-known general powerful method to synthesis of modular systems in various 
domains (e.g., [3], [26], [33], [66], [77]). MA is based on divide and conquer technique. A hierarchical 
structure of the designed system is a basis for usage of the method. The following basic partitioning 
techniques can be used to obtain the required system hierarchical model: (a) partitioning by system 
component/parts, (b) partitioning by system functions, (c) partitioning by system properties/attributes, 
and (d) integrated techniques. In this article, system hierarchy of system components (parts, subsystems) 
is considered as a basic one. This case corresponds to modular systems which are widely used in many 
domains of engineering, information technology, and management (e.g., [4], [24], [32], [33], [36], [62] . 
[75]). Many years the usage of morphological analysis in system design was very limited by the reason 
that the method leads to a very large combinatorial domain of possible solutions. On the other hand, 
contemporary computer systems can solve very complex computational problems and hierarchical system 
models can be used as a basis for partitioning/decomposition solving frameworks. Recent trends in the 
study, usage, and modification/extension of morphological analysis may be considered as the following: 
(1) hierarchical systems modeling, (2) optimization models, (3) multicriteria decision making, and (4) 
taking into account uncertainty (i.e., probabilistic and/or fuzzy estimates). 

In the article, the following system design methods are briefly described: (i) basic morphological anal- 
ysis (as morphological generation of admissible composite solutions) , (ii) modification of MA as method 
of closeness to ideal point (s), (iii) reducing of morphological analysis to optimization model as linear 
programming, (iv) multiple choice problem, (v) quadratic assignment problem, (vi) multicriteria analysis 
of morphological decisions with revelation of Pareto-efficient solutions, (vii) Hierarchical Morphological 
Multicriteria Design (HMMD) approach, and (viii) version of Hierarchical Morphological Multicriteria 
Design (HMMD) approach based on fuzzy estimates. The above-mentioned methods are illustrated by 
solving schemes, mathematical models, and illustrative numerical examples. Preliminary materials (a 
description of the morphological methods and an example for GSM network) were published in |41j , |50] . 
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2. Basic Configuration Problem 

Generally, morphological system design approaches are targeted to design of system configuration 
as a selection of alternatives for systems parts (e.g., [39]). Fig. 1 illustrates this problem. Here a 
composite (modular) system consists of m system parts: {-P(l), P(i), • P(m)}. For each system part 
(i.e., Vi, i — l,m) there are corresponding alternatives {X\ 1 X^, X l qi }, where qi is the number of 
alternatives for part i. Thus, the problem is: 

Select an alternative for each system part while taking into account some local and/or global objec- 
tives/preferences and constraints. 

Evidently, the objective/prereferences and constraints are based on (correspond to) quality of the 
selected alternatives and quality of compatibility among the selected alternatives. In [35], some other 
system configuration problems are described as well (e.g., reconfiguration, selection and allocation). 



System parts: { P(l), P(i), P(m) } 
System configuration example: 

S\ = X*2 "A- ... * A| -A- ... * X™ 
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Fig. 1. System configuration problem (selection) 

3. Morphological Design Approaches 

Our basic list of morphological design approaches consists of the following: (1) the basic version of 
morphological analysis (MA) (e.g., [7J, PH], [25], [SI]) [77] ) ; (2) the modification of morphological analysis 
as searching for an admissible (by compatibility) element combination (one representative from each mor- 
phological class, i.e., a set of alternatives for system part/component) that is the closest to a combination 
consisting of the best elements (at each morphological class) (e.g., [3J, [13]); (3) modification of morpho- 
logical analysis via reducing to linear programming (MA&linear programming) }29j : (4) modification of 
morphological analysis via reducing to multiple choice problem (MCP) or multicriteria multiple choice 
problem (e.g., [55], [75]); (5) modification of morphological analysis via reducing to quadratic assignment 
problem (QAP) (e.g., [11], [33], [39]); (6) the multicriteria modification of morphological analysis as fol- 
lows (Pareto-based MA): (a) searching for all admissible (by compatibility) elements combinations (one 
representative from each morphological class), (b) evaluation of the found combinations upon a set of 
criteria, and (c) selection of the Pareto-efficient solutions (e.g., [14], [18]); (7) hierarchical morphological 
multicriteria design (HMMD) approach ([33], [55]): and (8) version of hierarchical morphological multi- 
criteria design approach with probabilistic and/or fuzzy estimates (HMMD&uncertainty) [33]. Table 1 
contains some properties of the approaches above. 



Table 1. Description of approaches 



Method 


Scale for DAs 


Scale for IC 


Quality for decision 


References 


l.MA 


None 


{0,1} 


Admissibility 


[25], [55], [77] 


2. Ideal-point method 


None 


{0,1} 


Closeness to ideal point 




3. MA&linear programming 


Quantitative 


{0,1} 


Additive function 


m 


4. Multiple choice problem 


Quantitative 


None 


Additive function 


[56], [70] 


5. QAP 


Quantitative 


Quantitative 


Additive function 




6. Pareto-based MA 


None 


{0,1} 


Multicriteria 


[m, m 


7. HMMD 


Ordinal 


Ordinal 


Multicriteria 


[33], [36] 


8.HMMD & uncertainty 


Ordinal/fuzzy 


Ordinal/fuzzy 


Multicriteria 


[33] 
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3.1. Morphological Analysis 

The MA approach consists of the following stages: 

Stage 1. Building a system structure as a set of system parts/components. 

Stage 2. Generation of design alternatives (DAs) for each system part (i.e., a morphological class). 

Stage 3. Binary assessment of compatibility for each DAs pair (one DA from one morphological class, 
other DA from another morphological class). Value of compatibility 1 corresponds to compatibility of 
two corresponding DAs, value corresponds to incompatibility. 

Stage 4- Generation of all admissible compositions (one DA for each system part) while taking into 
account compatibility for each two DAs in the composition obtained. 

The method above is an enumerative one. Fig. 2 illustrates MA (binary compatibility estimates are 
depicted in Table 2). 



Table 2. Binary compatibility 




Fig. 2. Illustration for MA 
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Here the following morphological classes are examined: (a) morphological class 1 : { X{ , X\ , A3 , X\ , X\ } , 
(b) morphological class i\ {X{, X % 2 , A3, A|, A5}, and (c) morphological class m: {XT, XT, XT}- Further 
a simplified case is considered for three system parts (and corresponding morphological classes). The re- 
sultant (admissible) solution (composition or composite design alternative) is: Si = X\ *...*A|*...*A™. 

3.2. Method of Closeness to Ideal Point 

First, modification of MA as method of closeness to ideal point was suggested (e.g., [3], [13]). Illustration 
for method of closeness to ideal point is shown in Fig. 3 (binary compatibility estimates are contained 
in Table 3). 



Table 3. Binary compatibility 



Examples: S\ - 


= xl * . 
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Fig. 3. Illustration for MA with ideal point 
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Here for each system part (from the corresponding morphological class) the best design alternatives 
(as an ideal) are selected (e.g., by expert judgment). In the illustrative example (Fig. 3), the ideal design 
alternatives are: X\, A3, and A™. Thus, the ideal point (i.e., solution) is: So — X\ * ... ★ X\ ★ ... ★ A™. 
Unfortunately, this solution S is inadmissible (by compatibility). Admissible solutions are the following: 
Si = A 2 1 *...*A^*...*A™ and S 2 = Ag 1 * ... * X\ * ... * A™. 
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Let p(S',S") be a proximity (e.g., by elements) for two composite design alternatives S',S" G 
{S}. Then it is reasonable to search for the following solution S* G {S a } C {S} ({S a } is a set 
of admissible solutions): S* — Arg minjgK.i p(S,S Q ). Clearly, in the illustrative example solution 
£2 = Xl * ... *X^ * ... *X™ is more close to ideal solution S (i.e., p(S<2, S ) ^ p{S%, S ))- Generally, 
various versions of proximity (as real functions, vectors, etc.) can by examined (e.g., [3J, |13|). 

3.3. Pareto-Based Morphological Approach 

An integrated method (MA and multicriteria decision making, an enumerative method) was suggested 
as follows (e.g., [4]): 

Stage 1. Usage of basic MA to get a set of admissible compositions. 

Stage 2. Generation of criteria for evaluation of the admissible compositions. 

Stage 3. Evaluation of admissible compositions upon criteria and selection of Pareto-efficient solutions. 

Fig. 4 illustrates Pareto-based MA. Concurrently, binary compatibility estimates are depicted in Table 
4. Here admissible solutions are the following: Si = X^-k-.-kX^-k.-.-kXi 1 , S2 = X\-k ...+XI* ...*X™ , and 
S3 = Xl * ... * XI * ... *X™. Further, the solutions have to be evaluated upon criteria and Pareto-efficient 
solution(s) will be selected. It is important to note that the estimate vector for each DA can contain 
estimates of compatibility as well. Pareto-based morphological approach was used by several students 
during the author's course (instead of HMMD) ([37], gO], 03J). 



Table 4. Binary compatibility 
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Fig. 4. Illustration for Pareto-based MA 
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3.4. Linear Programming 

In [29j . MA is reduced to linear programming. Here constraints imposed on the solution are reduced 
to a set of inequalities of Boolean variables and quality criterion for the solution as an additive function 
is used. In this case, well-known methods for linear programming problems can be used. 

3.5. Multiple Choice Problem 

The basic knapsack problem is (e.g., [TH], [37J, [SO] - ): 

m m 

max CjXj s.t. < 6, Xi G {0,1}, i = l,m, 

i=l i=l 

where Xi = 1 if item i is selected, Cj is a value ("utility") for item i, and a, is a weight (or resource 
required). Often nonnegative coefficients are assumed. The problem is NP-hard ([H], [5D]) and can be 
solved by enumerative methods (e.g., Branch-and-Bound, dynamic programming), approximate schemes 
with a limited relative error (e.g., [37J, [SO])- In multiple choice problem (e.g., [50]), the items are divided 
into groups and we select element (s) from each group while taking into account a total resource constraint 
(or constraints). Here each element has two indices: where i corresponds to number of group and 

j corresponds to number of item in the group. In the case of multicriteria description, each element (i.e., 
has vector profit: Ci~j — (c\ j, cfj, c\J) and multicriteria multiple choice problem is |56| : 

m qi m qi qi 

max c i-j x ij 1 V£= l,r s.t. ajjXjj < b, Xjj = 1 Vi = l,m, iy G {0,1}. 

»=i j=i »=i j=i i=i 
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For this problem formulation, it is reasonable to search for Parcto-efBcient solutions. This design approach 
was used for design and redesign/improvement of applied systems (software, hardware, communication) 
( |44) . |56j . [70]). Here the following solving schemes can be used [SB]: (i) enumerative algorithm based 
on dynamic programming, (ii) heuristic based on preliminary multicriteria ranking of elements to get 
their priorities and step- by-step packing the knapsack (i.e., greedy approach), (iii) multicriteria ranking 
of elements to get their ordinal priorities and usage of approximate solving scheme (as for knapsack) 
based on discrete space of system excellence (as later in HMMD). 

3.6. Quadratic Assignment Problem 

Assignment/allocation problems are widely used in many domains (e.g., |19j . |63j). Simple assign- 
ment problem involves nonnegative correspondence matrix T = ||cij| (i = 1, n, j = 1, n) where Cij is 
a profit ('utility') to assign element i to position j. The problem is (e.g., |19j): 

Find assignment n = (7r(l), 7r(n)) of elements i (i = 1, n) to positions ir(i) which corresponds 

to a total effectiveness: Ci j7r (;) — > max. 

More complicated well-known model as quadratic assignment problem (QAP) includes interconnection 
between elements of different groups (each group corresponds to a certain position) (e.g., [llj . [33], [63]). 
Let a nonnegative value d(i, j\ , fc, j'2) be a profit of compatibility between item j\ in group Ji and item j'2 
in group Jk- Also, this value of compatibility is added to the objective function. QAP may be considered 
as a version of MA. Thus, QAP can be formulated as follows: 

m qt qi qu 

max ^2^2 Ci d x hj +^2 X! X! d ( l '3i-' k 'h) x l,ji x k,fr, 1 = 1, m, k=l ) m; 

1=1 j=i i < kj 1 =ij 2 =i 

m qi qi 

s ~t- Xi >i — ^' ^2 Xi -j -5 1 V£ = 1, m, Xi.j £ {0, 1}. 

i=i j=i j=i 

QAP is NP-hard. Enumerative methods (e.g., branch-and-bound) or heuristics (e.g., greedy algorithms, 
tabu search, genetic algorithms) are usually used for the problem. In the case of multicriteria as- 
signment problem, the objective function is transformed into a vector function, i.e., Cjj =>■ cTJ = 
(cj j, cf ■, c\ j) and the vector objective function is, for example: 

m n m n m n 

i—l j—1 i—1 j—1 i—1 j—1 

Here Pareto-efficient solutions are usually searched for. 

3.7. Hierarchical Morphological Multicriteria Design 

A basic description of Hierarchical Morphological Multicriteria Design (HMMD) is contained in ( |33] . 
[36] , [39]). The assumptions of HMMD are the following: (a) a tree- like structure of the system; (b) 
a composite estimate for system quality that integrates components (subsystems, parts) qualities and 
qualities of interconnections (hereinafter referred as 'IC') across subsystems; (c) monotonic criteria for 
the system and its components; and (d) quality of system components and IC are evaluated on the basis 
of coordinated ordinal scales. The designations are: (1) design alternatives (DAs) for nodes of the model; 
(2) priorities of DAs (r — 1, fc; 1 corresponds to the best level); (3) ordinal compatibility estimates for 
each pair of DAs (w = 0,1; I corresponds to the best level). The basic phases of HMMD are: 

Phase 1. Design of the tree-like system model (a preliminary phase). 

Phase 2. Generating DAs for leaf nodes of the system model. 

Phase 3. Hierarchical selection and composing of DAs into composite DAs for the corresponding higher 
level of the system hierarchy (morphological clique problem) . 

Phase 4- Analysis and improvement of the resultant composite DAs (decisions). 

Let S be a system consisting of m parts (components): P(l), P(i), P(m) (Fig. 1). A set of 
DAs is generated for each system part above. The problem is: 

Find composite design alternative S = S(l) * S(i) S(m) of DAs (one representative design 
alternative S(i) for each system component/part P(i), i = 1, m) with non- zero IC estimates between 
the selected design alternatives. 
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A discrete space of the system quality is based on the following vector (Fig. 5): N(S) = (w(S); n(S)), 
where w(S) is the minimum of pairwise compatibility between DAs which correspond to different system 
components (i.e., V Pj 1 and Pj 2 , 1 < j% ^ ji < m) in S, n(S) — (ni, n r , ■■•Tik), where n r is the 
number of DAs of the rth quality in S (X) r =i T7-r = TO )' Here composite solutions (composite DAs) are 
searched for which are nondominated by N(S) (i.e., Pareto-efHcient solutions) (Fig. 5). 

In ([33], [35]), the described combinatorial problem is called morphological clique problem, this problem 
is NP-hard (because a more simple its subproblem is NP hard [2S]). Generally, the following layers of 
system excellence can be considered (e.g., [33] ): (i) ideal point; (ii) Pareto-efHcient points; (hi) a 
neighborhood of Pareto-efHcient DAs (e.g., a composite decision of this set can be transformed into a 
Pareto-efHcient point on the basis of an improvement action(s)). Clearly, the compatibility component of 
vector N(S) can be considered on the basis of a poset-like scale too (as n(S) ). In this case, the discrete 
space of system excellence will be an analogical lattice ([31], |36|). 

Fig. 6 and Fig. 7 illustrate HMMD (by a numerical example for a system consisting of three parts 
S = X *Y -k Z). Priorities of DAs are shown in Fig. 6 in parentheses and are depicted in Fig. 7. Table 5 
contains compatibility estimates (they are pointed out in Fig. 7 too). In the example, composite decisions 
are (Pareto-efHcient solutions) (Fig. 5, Fig. 6, Fig. 7, Fig. 8): Si = X 2 *Y 2 *Z 2 , N(Si) = (1;2,1,0); 
S 2 =Xi*Y 2 *Z 2 , JV(5 2 ) = (2;1,2,0); S 3 = X 1 *Y 1 *Z 3 , N(S 3 ) = (3; 0, 2, 1). 

HMMD method was used to design various modular systems (e.g., packaged software, communication 
networks, security system, web-hosting system, car, telemetric system, test inputs in system testing, 
concrete technology, immunoassay technology, management for smart home, sensor node) ( 33 , [35 , 36 , 
EE], 05, 05J, 0BJ, g5], 05], [51], [5J, [S5J, [SB]). In addition, HMMD was used in modular system 
improvement processes ([55 ] , [5S ] . [5B]. 02 ] . 07]. [57] V 
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Fig. 5. Space of system quality (3 system parts, 3 levels of element quality) 
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Fig. 6. Example of composition Fig. 7. Concentric presentation 




Table 5. Compatibility 
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Fig. 8. Illustration for space of quality 



3.8. Morphological Multicriteria Design (Uncertainty) 

The version of HMMD under uncertainty has been suggested in [33]. Here a brief description of the 
approach is presented (by a simplified example). Fuzzy estimates for DAs or/and IC are considered. The 
following designations are used: (1) i is an index corresponding to the design alternative; (2) /i r (i) is 
a membership function of priority r(i) G {1,2,3}, we consider the following set: fJq(i), I — 1,...,3 (I 
corresponds to value of priority); and (3) fi w (ix,i2) is a membership function of compatibility w(ix,i2) 
between alternatives i\ and Z2, we use the following set: {^fc(*i, 22), k = 0, ...,3} (k corresponds to 
value of pairwise compatibility). Thus, let fj, w (i±,i2) be the following vector (i.e., pairwise compatibility 
vector): i 2 ), (J% {kM), Vi(h, h), Mo ^))- 

Now let r a (i) and w a (11,12) be aggregated estimates for design alternative i, and for a pair of design 
alternatives (11,12) accordingly. 

Here a basic system example is depicted in Fig. 9. Table 6 contains normalized fuzzy priorities {fJ-l(i)} 
and realistic aggregated priorities {r a (i)}. Table 7 contain normalized fuzzy compatibility {/i^ 12)}, 
and Table 8 presents realistic aggregated compatibility {ty Q («i, 22)}- 
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Table 6. Fuzzy priorities 



DAs i 




M2OO 


M3OO 


r a (i) 


Ax 


1.00 


0.00 


0.00 


1 


A 2 


0.00 


0.05 


0.95 
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Bx 


0.15 


0.65 


0.20 


2 


B 2 


0.85 


0.15 


0.00 


1 


Cx 


1.00 


0.00 


0.00 


1 


c 2 


0.00 


1.00 


0.00 


2 



Thus, the following relationships over the set of DAs can be obtained (by fuzzy priorities): (a) Ax > A2; 
(b) Bx > B2; and (c) C\ and C2 are non-comparable. Now it is reasonable to examine the following 
cases: 

Case 1: deterministic (aggregated) estimates of priorities for DAs {r a (i)} and deterministic (aggre- 
gated) estimates of compatibility for IC {u> Q (ii,i2)} (a basic case). 

Case 2: estimates of DAs are aggregated (deterministic) {r a (i)}, and estimates of compatibility are 
fuzzy {^(ii,i 2 )}, V(ii,i 2 ). 

Case 3: estimates of DAs are fuzzy {/i[(i)} Vi and estimates of IC are aggregated (deterministic) 
{w a (ii,i 2 )}, V(ii,i 2 ). 

Case 4- estimates of DAs are fuzzy Vi and estimates of compatibility are fuzzy {/ij? (ix, 12)}, 

V(ii,» 2 ). ' * " 
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Table 7. Fuzzy compatibility 
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0.5; 0.2; 0.3; 0.0 0.0; 0.3; 0.4; 0.3 0.0; 0.4; 0.5; 0.1 0.2; 0.4; 0.3; 0.1 
0.1; 0.2; 0.4; 0.3 0.7; 0.3; 0.0; 0.0 0.4; 0.4; 0.2; 0.0 0.0; 0.7; 0.3; 0.0 

0.6; 0.3; 0.1; 0.0 0.0; 0.5; 0.5; 0.0 
0.0; 0.4; 0.5; 0.1 0.2; 0.5; 0.3; 0.0 



Table 8. Aggregated compatibility 
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Fig. 10 illustrates the cases above (top index of composite DAs corresponds to the case). 

Clearly, that main solving method is based on two stages: (1) generation of feasible composite decisions; 
and (2) selection of Pareto-efficient decisions. Unfortunately, it is reasonable to point out the following 
two significant features of our synthesis problem with fuzzy estimates: (a) complexity of corresponding 
combinatorial problems is increasing because a number of analyzed composite decisions is more than in 
deterministic case; and (b) it is necessary to construct a preference rule to select the best fuzzy decision(s). 



Criterion 1: n(S) 
(quality by elements) 



The .ideal 
point 

Q 




Criterion 2: w(S) 
(quality by IC) 

Fig. 10. Illustrative space for composite DAs (from [33] ) 



4. Design Examples for GSM Network 

In recent two decades, the significance of GSM network has been increased (e.g., [12], [20], [22], [58], 
[6"T] . [64] . [74]). Thus, there exists a need of the design and maintenance of this kind of communication 
systems. Here a numerical example for design of GSM network (a modification of an example from [50] ) 
is used to illustrate and to compare several MA-based methods: basic MA, method of closeness to ideal 
point, Pareto-based MA, multiple choice problem, and HMMD. 
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4.1. Initial Example 

The general tree-like simplified model of GSM network is as follows (Fig. 11, the developers of DAs 
are pointed out in parentheses): 



0. GSM network S = A * B. 

1. Switching Subsystem SSS ( A = M * L), 

1.1. Mobile Switching Center/ Visitors Location Register MSC/VLR M : M x (Motorola), M 2 (Alca- 
tel), M 3 (Huawei), M4 (Siemens), and Ms (Ericsson). 

1.2. Home Location Register/Authentification Center HLR/AC L : L\ (Motorola), L 2 (Ericsson), L 3 
(Alcatel), and L 4 (Huawei). 

2. Base Station SubSystem BSS (B = V * U * T). 

2.1. Base Station Controller BSC V : V x (Motorola), V 2 (Ericsson), V 3 (Alcatel), V4 (Huawei), V 5 
(Nokia), and Vq (Siemens). 

2.2. Base Transceiver Station BTS U : U\ (Motorola), U 2 (Ericsson), U 3 (Alcatel), U4 (Huawei), and 
U s (Nokia). 

2.3. Transceivers TRx T : T t (Alcatel), T 2 (Ericsson), T 3 (Motorola), T 4 (Huawei), and T 5 (Siemens). 
Note an initial set of possible composite decisions contained 3000 combinations (5x4x6x5x5). 



GSM network S = A * B = (M * L) * (V * U * T) 



SSS A = M*L 



MSC/VLR M: M x , 
M 2 ,M 3 ,M A ,M 5 



HLR/AC L: 
L\,L 2l L 3l L4 



BSS B = V*U*T 










BSS V: Vi,V 2 , 

v 3 ,v 4 ,v 5 ,v 6 


BTS U: U lr 

U 2 ,U 3 ,U 4 ,U 5 


TRx T: Ti, 
T 2 ,T 3 ,T 4 ,T 5 



Fig. 11. General simplified structure of GSM network 

Let us consider criteria for system components as follows (weights of criteria are pointed out in paren- 
theses): 

1. M: maximal number of datapathes (1000 pathes) (C m i, 0.2); maximal capacity VLR (100000 
subscribers) (C m2 , 0.2); price index ( 100000/price(USD) ) (C m3 , 0.2); power consumption ( 1/power 
consumption (kWt) ) (C m4l 0.2); and number of communication and signaling interfaces (C m 5, 0.2). 

2. L: maximal number of subscribers (100000 subscribers) (C;i, 0.25); volume of service provided 
(C /2 , 0.25); reliability (scale [1, 10]) (C /3 , 0.25); and integratability (scale [1, 10]) (C M , 0.25 ). 

3. V: price index ( 100000/cost (USD) ) (C„i, 0.25); maximal number of BTS (C v2 , 0.25); handover 
quality (C v3 , 0.25); and throughput (C v4 , 0.25). 

4. U: maximal number of TRx (C uU 0.25); capacity (C u2 , 0.25); price index ( 100000/cost(USD) ) 
(C u3 , 0.25); and reliability (scale [1, 10]) (C„ 4 , 0.25). 

5. T: maximum power-carrying capacity (Cti, 0.3); throughput (Ct 2 , 0.2); price index ( 100000/ 
cost(USD) ) (C t3 , 0.25); and reliability (scale [1, 10]) (C t4 , 0.25). 

Tables 9, 10, 11, 12, and 13 contain estimates of DAs upon criteria above (data from catalogues, expert 
judgment) and their resultant priorities (the priorities are based on multicriteria ranking by an outranking 
technique [67J. Compatibility estimates are contained in Tables 14 and 15 (expert judgment). 

Table 9. Estimates for M 



DAs 


Criteria 


Priority 
r 






C m3 




Cm 5 


Mi 


3.7 


8.6 


6 


5.1 


4 


2 


M 2 


4.0 


11 


8 


7 


5 


3 


M 3 


4.1 


10 


9 


7 


4 


3 


M A 


3.2 


7 


5 


6 


3 


1 


M s 


3.5 


8.7 


6.2 


5 


4 


2 



10 



Table 11. Estimates for V 



Table 10. Estimates for L 



DAs 


Criteria 


Priority 




C n 








r 


Li 


9 


7 


7 


8 


1 


L 2 


10 


4 


9 


8 


1 


L 3 


12 


8 


10 


10 


2 


Li 


9 


5 


8 


8 


1 


Table 12. Estimates for U 


DAs 


Criteria 


Priority 






C u2 


C U 3 


C u4 


r 




2 


7 


5 


8 


1 


u 2 


4 


10 


6 


10 


3 


u 3 


3 


9 


6 


10 


2 


Ui 


3 


6 


3 


7 


1 


u 5 


3 


10 


6 


9 


2 



Table 14. Compatibility 







L 2 


L 3 


u 


Mi 


3 


2 





3 


M 2 


2 


3 


2 





M 3 





2 


3 


2 


M 4 


2 


3 


3 


3 


M 5 


3 


3 





3 



DAs 


Criteria 


Priority 




C v i 


C V 2 


C v3 




r 


Vi 


6 


4 


3 


4 


1 


v 2 


7 


5 


7 


7 


2 


v 3 


9 


7 


10 


7 


3 


v 4 


7 


5 


8 


6 


2 


v 5 


6 


3 


4 


4 


1 


v e 


10 


6 


9 


7 


3 


Table 13. Estimates for T 


DAs 


Criteria 


Priority 




C tl 


Cfi 


c t3 




r 


Ti 


9 


7 


10 


7 


3 


T 2 


6 


4 


3 


4 


1 


T 3 


7 


5 


7 


7 


2 


T 4 


7 


5 


8 


6 


2 


T 5 


6 


3 


4 


4 


1 



Table 15. Compatibility 







u 2 


u 3 


Ui 


c/ 5 


Ti 


T 2 


r 3 


Ti 


r 5 


Vi 


2 


2 


2 


2 


3 


3 


2 


2 


2 


2 


v 2 


3 


3 


3 


2 








3 





3 


2 


v 3 


3 


3 


3 


2 








3 





3 


2 


V 4 


3 


2 





2 


3 





2 





2 


2 


v 5 


3 








2 





2 


2 





2 


2 


v e 





3 


2 


3 


2 


3 





2 


2 

















2 








2 


3 


u 2 















2 





3 





u 3 















2 





3 





Ui 















3 


3 








u 5 












3 





2 


2 






4.2. Morphological Analysis 

In the case of basic MA, binary compatibility estimates are used. To decrease the dimension of the 
considered numerical example, the following version of MA is examined. Let us consider more strong 
requirements to compatibility: (Tables 16 and 17): (i) new compatibility estimate equals 1 if the old 
estimate was equal 3, (ii) new compatibility estimate equals 1 if the old estimate was equal or 1 or 
2. Clearly, here we can get some negative results, for example: (a) admissible solutions are absent, 
(b) some sufficiently good solutions (e.g., solutions with one/two compatibility estimate at the only 
admissible/good levels as 1 or 2) will be lost. As a result, the following admissible DAs can be analyzed: 

(1) nine DAs for A: A 1 = M x * Li, A 2 = Mi * L 4 , A 3 = M 2 * L 2 , A 4 = M 3 * L 3 , A 5 = M 4 * L 2 , 
A 6 = M 4 * L 3 , A 7 = M 5 * Li, A 8 = M 5 ★ L 2 , and A g = M 5 * L 4 ; 

(2) five DAs for B: B^. =V 1 *U b *T l , B 2 = V 2 *U 2 * T 4 , B 3 = V 2 *U 3 * T 4 , B 4 = V 3 *U 2 * T 4 , and 
B 5 = V 3 *U 3 *Ti; 

and the resultant composite DAs are: Si — A\ * Bi, S 2 
S5 = A§ * Bi , Sq — Aq -k Bi , 57 = A-; -k B x , 5s = As * B x , So, 

Sn — A§-kB 2 , S15 = Aq*B 2 , Sie 
S 2 i = A 3 -kB 3 , S 22 = Ai*B 3 , S 23 = As-kB 3 , 
S 2 s = Ai-kBi, S 2 g = A 2 -kBi, S 3 q = A 3 -kBi, 

5*35 = A 8 *Bi 7 S 3 q = Ag-kBi, S 3 j = Ai-kBs, 
S42 = Aq * B$, 5*43 = A7 -k i? 5 , S44 = As* _B 5 , 



Si 2 

S19 

S 2 a 

s 33 

Sio 



A 3 *B 2 , 
Ai-kB 3 , 
A S *B 3 , 
A 6 *Bi, 
Ai*B 5: 



S13 = Ai-kB 2 , 
S 2 o = A 2 *B 3 , 
S 2 7 = Ag -k B 3 ; 
S 3 i = Aj -k Bi, 
Sn = A$ -k B5 , 



= A 2 kB u S 3 = A 3 *B U 
Ag ★ Bi\ Sio = Ai ★ B 2 , 
A 7 *B 2 , S17 = A S *B 2 
S 2i = A 6 -kB 3 
S 3 i = AikBi 

S 3 8 = A 2 kB 5 

and 545 = A 



Si — Ai -k Bi , 
5*11 = A 2 k B 2 , 
Sis — Ag-kB 2 ; 
S 25 =A 7 kB 3 , 
S 32 = A 5 kBi, 

S 3 g = A 3 kB 5 , 

kB 5 . 
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Finally, the next step has to consist in selection of the best solution. 

Table 17. Compatibility 



Table 16. Compatibility 







L 2 


L 3 


L 4 


Mi 


1 








1 


M 2 





1 








M 3 








1 





M A 





1 


1 


1 


M 5 


1 


1 





1 







u 2 


u 3 


u 4 


u 5 


Ti 


T 2 


T 3 


T 4 


T 5 


Vi 














1 


1 














v 2 


1 


1 


1 











1 





1 





v 3 


1 


1 


1 











1 





1 





v 4 


1 











1 

















v 5 


1 





























V e 





1 





1 





1 






































1 


u 2 





















1 





u 3 





















1 





Ui 















1 


1 








u 5 












1 















4.3. Method of Closeness to Ideal Point 

Here the initial set of admissible solutions corresponds to the solution set that was obtained in previous 
case (i.e., basic MA). Evidently, this approach depended on the kind of the proximity between the ideal 
point ( S 1 ) and examined solutions. 

First of all, let us consider estimate vector for each admissible solution (basic estimates are contained 
in Tables 9, 10, 11, 12, and 13): 

z = (z M [J z L [J z v [J zu [J z T ) = 

(z m i,z m2 , Z m3 , Zmi, Z m 5i Z[i, Zl 2 , Zi 3 , Zl 4 , Z v \, 

^v2i ^v3i %v4-> Z>ul-> %u2i ^u3i Zui , Ztl , Zfl , Zt 3 , Ztij- 

On the other hand, it may be reasonable to consider a simplified version of the estimate vector as follows: 
z = (rj\/, r\L, ry , ru, tt), where tm , f\L, fV ,ru ,Tt are the priorities of DAs which are obtained for local 
DAs (for M, for L, for V, for U, and for T; Table 9, Table 10, Table 11, Table 12, Table 13). To simplify 
the considered example, the second case of the estimate vector is used. Thus, the resultant vector 
estimates (i.e., {z} ) for examined 45 admissible solutions are contained in Tables 18 and 19. Evidently, 
it is reasonable to consider the estimate vector for the ideal solution as follows: z} — (1,1,1,1,1). Now 
let us use a simplified proximity function between ideal solution / and design alternative DA as follows 
(i.e, metric like I 2 ): 



p(I,DA)= Yl (z k (I)-z k (DA)f. 

Y ke{M,L,V,U,T} 

The resultant proximity is presented in Tables 18 and 19. Finally, the best composite design alternative 
(by the minimal proximity) is: 5^ = 523 = A$ * B 3 = M 3 * L\ * Vi * U 2 * T 3 (p = 1.7321). Several 
composite DAs are very close to the best one, for example: 

S( = S w = A!*B 3 = M^L^V^U^Ti (p = 2.0), S% = S 2Q = A 2 *B 3 = M^L^V^U^Ti (p = 2.0), 
Si - S 24 - A e *B 3 = M 4 *L 3 *V 2 *U 3 *T 4 (p = 2.0), Sj = S 25 = A 7 *B 3 = M s * L^V 2 *U 3 ^T 4 (p = 2.0), 
Si = S 2e = A S *B 3 = M 5 * L 2 * V 3 * U 2 * T 4 (p = 2.0) , and S£ = S 27 = A 9 *B 3 = M 5 * L i *V 2 *U s *T 4 
(p = 2.0). 

It may be reasonable to point out several prospective directions for the improvement of this method: 

(1) consideration of special types of proximity between solutions and the ideal point (e.g., ordinal 
proximity, vector-like proximity [33] , etc.); 

(2) usage of special expert judgment interactive procedures for the assessment of the proximity; 

(3) consideration of a set of ideal points (the set can be generated by domain expert (s)); and 

(4) design of special support visualization tools which will aid domain expert (s) in his/her (their) 
activity (i.e., generation of the ideal point and assessment of proximity). 

In addition let us list the basic approaches to generation of the ideal point (s): 
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1 . consideration of design alternative with the estimate vector in which each component equals the best 
value of the design alternatives estimates (by the corresponding criterion, i.e., minimum or maximum); 

2. consideration of design alternative with the estimate vector in which each component equals the 
best value of the corresponding criterion scale (i.e., minimum or maximum); 

3. expert judgment based generation design alternative(s); 

4. projection of expert judgment based design alternatives into convex shell of the set of Pareto-efficient 
points; etc. 



Table 18. Estimates of admissible solutions 



DAs 


z 


Closeness to 
ideal point 


Membership 
of Pareto-set 


Si 




1 
1 


1 

1 


z 


O ) 




Mr, 
IN O 


So 


(2, 


1 


1 


2 


3) 




1 N O 


s 3 


(3, 


1, 


1 


2 


3) 


3.0 


No 


s 4 


(3, 


2 


1 


2 


3) 


3.1623 


No 


s 5 


(1, 


1 


1 


2 


3) 


2.2361 


Yes 


s 6 


(1, 


2 


1 


2 


3) 


2.4495 


No 


s 7 


(2, 


1, 


1 


2 


3) 


2.4495 


No 


s$ 


(2, 


1 


1 


2 


3) 


2.4495 


No 


s 9 


(2, 


1 


1 


2 


3) 


2.4495 


No 


Sw 


(2, 


1 


2 


3 


2) 


2.6458 


No 


Su 


(2, 


1, 


2 


3 


2) 


2.6458 


No 


Sl2 


(3, 


1 


2 


3 


2) 


3.1623 


No 


Sl3 


(3, 


2 


2 


3 


2) 


3.3166 


No 


5*14 


(1, 


1 


2 


3 


2) 


2.4495 


No 


5*15 


(1, 


2 


2 


3 


2) 


2.6458 


No 


Si 6 


(2, 


1 


2 


3 


2) 


2.6458 


No 


S17 


(2, 


1 


2 


3 


2) 


2.6458 


No 


Sis 


(2, 


1, 


2 


3 


2) 


2.6458 


No 


S19 


(2, 


1, 


2 


2 


2) 


2.0 


No 


S20 


(2, 


1 


2 


2 


2) 


2.0 


No 


S21 


(3, 


1 


2 


2 


2) 


2.6458 


No 


S22 


(3, 


2 


2 


2 


2) 


2.8284 


No 


S23 


(1, 


1, 


2 


2 


2) 


1.7321 


Yes 



Table 19. Estimates of admissible solutions 



DAs 


z 


Closeness to 
ideal point 


Membership 
of Pareto-set 


S24 


(1 


2 


2 


2,2) 


2.0 


No 


S25 


(2 


1 


2 


2,2) 


2.0 


No 


S26 


(2 


1 


2 


2,2) 


2.0 


No 


S27 


(2 


1 


2 


2,2) 


2.0 


No 


S28 


(2 


1 


3 


3,2) 


3.1623 


No 


S29 


(2 


1 


3 


3,2) 


3.1623 


No 


S30 


(3 


1 


3 


3,2) 


3.6056 


No 


S31 


(3 


2 


3 


3,2) 


3.7417 


No 


S32 


(1 


1 


3 


3,2) 


3.0 


No 


S33 


(1 


2 


3 


3,2) 


3.1623 


No 


S34 


(2 


1 


3 


3,2) 


3.1623 


No 


S35 


(2 


1 


3 


3,2) 


3.1623 


No 


S36 


(2 


1 


3 


3,2) 


3.1623 


No 


S37 


(2 


1 


3 


2,2) 


2.6458 


No 


S38 


(2 


1 


3 


2,2) 


2.6458 


No 


S39 


(3 


1 


3 


2,2) 


3.1623 


No 


S40 


(3 


2 


3 


2,2) 


3.3166 


No 


S41 


(1 


1 


3 


2,2) 


2.4495 


No 


S42 


(1 


2 


3 


2,2) 


2.6458 


No 


S43 


(2 


1 


3 


2,2) 


2.658 


No 


S44 


(2 


1 


3 


2,2) 


2.6458 


No 


S45 


(2 


1 


3 


2,2) 


2.6458 


No 



4.4. Pareto-based Morphological Analysis 

Here the initial set of admissible solutions corresponds to the previous case of basic MA. Two approaches 
can be used for mulricriteria assessment of admissible solutions: 

1. Basic method: selection of Pareto-efficient solutions over the set of admissible composite solutions 
on the basis of of usage of the initial set of criteria for assessment of each admissible composite DAs; 

2. Two-stage method: 

(i) assessment of initial components by the corresponding criteria and ranking of the alternative com- 
ponents to get an ordinal priority for each component, 

(ii) selection of Pareto-efficient solutions over the set of admissible composite solutions on the basis of 
of usage of the vector estimates which integrate priorities of solution components above. The results of 
the Pareto-based MA are presented in Tables 18 and 19, i.e., the resultant (Pareto-efficient) DAs are: (i) 
S[ = S 5 = A 5 -kB 1 = M i *L 2 *V 1 *U 5 -kT 1 and (ii) S% = S 23 = A 5 *B 3 = M 4 * L 2 *V 2 *U 3 *T 4 . 

4.5. Multiple Choice Problem 

Multiple choice problem with 5 groups of elements (i.e., for M, L, V, U, T) is examined (Fig. 12). 
Evidently, here it is reasonable to examine multicriteria multiple choice problem. In the example, a 
simplified problem solving approach is considered (Table 20): 

(i) a simple greedy heuristic based on element priorities is used; 

(ii) for each element (i.e., i, j) 'profit' is computed as follows: Cjj = 4 — nj; 

(iii) for each element (i.e., a required resource is computed as follows: a^j = 11 — Zij where Zij 
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equals: (a) for M: the estimate upon criterion C m3 (Table 9), (b) for L: 1.0, (c) for V: the estimate 
upon criterion C v \ (Table 11), (d) for U: the estimate upon criterion C u3 (Table 12), and (e) for T: the 
estimate upon criterion C mt3 (Table 13). 



• S = M*L*V *U *T 



M 


.1 


U 


iu 


It 


MSC/ 


HLR/ 


BSC 


BTS 


TRx 


VLR 


AC 


^i(l) 


Ui(l) 


Ti(3) 


Mi (2) 


Mi) 


y 2 (2) 


U 2 (S) 


r 2 (i) 


M 2 (3) 


L 2 (l) 


^3(3) 


Us(2) 


33(2) 


M 3 (3) 


i 3 (2) 


V 4 (2) 


U A {1) 


T 4 (2) 


M 4 (l) 


L 4 (l) 


^5(1) 


U 5 (2) 


r 6 (i) 


M 6 (2) 




^6(3) 







Fig. 12. Designed GSM network (priorities of DAs 
are shown in parentheses) 

Table 20. Example for multiple choice problem 



No. 


DAs 


Priority 


Resource 


c i,j / a i,j 


Selection 


Selection 






r 


requirement 

a i,j 




(constraint: 
< 14) 


(constraint: 
< 15) 


(1,1) 


Mi 


2 


5.0 


0.4 


No 


No 


(1,2) 


M 2 


3 


3.0 


0.33 


No 


No 


(1,3) 


M 3 


3 


2.0 


0.5 


No 


No 


(1,4) 


M 4 


1 


6.0 


0.5 


Yes 


Yes 


(1,5) 


M 5 


2 


4.8 


0.38 


No 


No 


(2,1) 


ii 


1 


1.0 


3.0 


Yes 


Yes 


(2,2) 


i 2 


1 


1.0 


3.0 


No 


No 


(2,3) 


L 3 


2 


1.0 


2.0 


No 


No 


(2,4) 


U 


1 


1.0 


3.0 


No 


No 


(3,1) 


Vi 


1 


5.0 


0.6 


No 


No 


(3,2) 


^2 


2 


4.0 


0.5 


No 


No 


(3,3) 


Vs 


3 


2.0 


0.5 


No 


No 


(3,4) 


^4 


2 


4.0 


0.5 


No 


No 


(3,5) 




1 


5.0 


0.6 


No 


No 


(3,6) 


^6 


3 


1.0 


1.0 


Yes 


Yes 


(4,1) 




1 


6.0 


0.5 


No 


Yes 


(4,2) 




3 


5.0 


0.2 


No 


No 


(4,3) 


u 3 


2 


5.0 


0.4 


Yes 


No 


(4,4) 


u 4 


3 


8.0 


0.39 


No 


No 


(4,5) 


u 5 


2 


5.0 


0.4 


No 


No 


(5,1) 


Ti 


3 


1.0 


1.0 


Yes 


Yes 


(5,2) 


T 2 


1 


8.0 


0.39 


No 


No 


(5,3) 


T 3 


2 


4.0 


0.5 


No 


No 


(5,4) 


T 4 


2 


3.0 


0.66 


No 


No 


(5,5) 


r 5 


1 


7.0 


0.42 


No 


No 



Thus, the following simplified one-objective problem is considered: 

5 qi 5 <2i qi 

maX X!X/ c ij x ij s -t- ^ ^ ClijXij < 6, Xij = 1 Mi = 1, 5, Xij G {0, 1}, 

i=l j=l i=l j = l j=l 

where qi = 5, q 2 = 4, q 3 = 6, g 4 = 5, q§ = 5. After the usage of the greedy heuristic, the following 
composite DAs are obtained (Table 20): 

(1) resource constraint b = 14: Sf = Af 4 * Li * V e * U 3 * T\, 

(2) resource constraint b = 15: = M 4 * Li * V e * U\ * T\. 
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4.6. Hierarchical Morphological Design 

A preliminary example for HMMD was presented in [50] (Fig. 13). For system part A, we get the 
following Pareto-efficient composite DAs: (1) A a = M 4 *L 2 , N(Ai) = (3; 2, 0, 0); (2) A 2 = M i -kL 4l , 
N(A2) = (3; 2, 0, 0). For system part £?, we get the following Pareto-efficient composite DAs: (1) B\ — 
V b *U x *T b , N(B Y ) = (2; 3, 0,0); (2) 5 2 = V 5 * f/ 4 * T 2 , N(B 2 ) = (2; 3, 0, 0); (3) fl 3 = ^ * £/- B *Ti, 
7V(s 3 ) = (3 ; 1,1,1), and (4) Bi = V 2 *U 3 * T 4 , 7V(5 4 ) = (3; 0, 3, 0). Fig. 14 illustrates system quality 
for B. 



Si - 

S 2 ~ 

S3 

5 4 : 

S 5 : 

S 6 : 

s 7 - 
sss 



= A-kB-- 
= A 1 *B 1 
= A X *B 2 
= At*B 3 
= A 2 *B 1 
= A 2 *B 2 
-A 2 *B 3 
= A 1 *B 4 
= A 2 *B 4 

A = M-k 



A! = 

A 2 = 



M 4 * L 2 
M4, * L 4 



= (M*L)*(V*U*T) 
= (Af 4 *L 2 )*(V 5 *£/i*T 5 ) 
= (M 4 *L 2 )*(V 5 *U 4 *T 2 ) 
= (M 4 *i2)*(Vi*l75*ri) 
= (Mi*Li)*(V 5 *Ui*T 5 ) 
= (M A * Li) * (V5 *Ui* T 2 ) 
= (Mi*Li)*(Vi*U 5 *T 1 ) 
= (Mi*L 2 )*(V 2 *U 3 *Ti) 
= (Mi*Li)*(V 2 *U 3 *Ti) 

BSS 

B! = 



L 



B =V*U*T 
V 5 *U!*T 5 
f B 2 = V 5 *U 4 *T 2 
B 3 = Vi*C/ 5 *Ti 
B 4 = V 2 *U 3 * T 4 



M 
MSC/ 
VLR 

Mr (2) 

M 2 (3) 
M 3 (3) 
M 4 (l) 
M 5 (2) 



71 

HLR/ 

AC 

ii(l) 
L 2 (l) 

M2) 
L 4 (l) 





\u 


It 


BSC 


BTS 


TRx 


^i(l) 


Ui(l) 


T x (3) 


y 2 (2) 


U 2 (3) 
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Fig. 13. Designed GSM network (priorities of DAs 
are shown in parentheses) 




Fig. 14. Space of system quality for B 



Now it is possible to combine the resultant composite DAs as follows (Fig. 13): 
(1) Sf = Ai *Bx = [Mi * L a ) * (Vfe *U X * T 5 ); (2) Sf = Ax * B 2 = (M 4 * L 2 ) * (V5 * t/ 4 * T 2 ); 
(3) = At * B 3 = (Mi * L 2 ) * (Vi *U 5 * Ti); (4) S i H = A 2 *B 1 = (M 4 * L 4 ) * (V5 * U x * T 5 ); 
(5) Sf = A 2 * B 2 = (M 4 * Li) * (V 5 *Ui* T 2 ); (6) S$ = A 2 * B 3 = (M 4 * L 4 ) * (Vi * J7 5 ★ Ti); 
(7) Sf = A x * B 3 = (Mi * L a ) * (Va * f/ 3 * and (8) Sf = A 2 * B 3 = (M 4 * L 4 ) * (V 2 *U 3 * T 4 ). 
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Finally, it is reasonable to integrate quality vectors for components A and B to obtain the following 
quality vectors: N{S?) = (2; 5, 0,0), N{S 2 H ) = (2; 5, 0,0), iV(Sf) = (3; 3, 1,1), AT(Sf) = (2; 5, 0,0), 
JV(Sf ) = (3; 3, 1,1), andiV(5f ) = (3; 3, 1,1). N(S?) = (3; 2, 3,0), and iV(5f ) = (3; 2, 3,0). Further the 
obtained eight resultant composite decisions can be analyzed to select the best decision (e.g., additional 
multicriteria analysis, expert judgment). 

4.7. Brief Comparison and Discussion of Methods 

Note, 45 resultant solutions were obtained by basic MA. Table 21 integrates resultant composite solu- 
tions for four methods: (1) ideal point method (the best solution and six close solutions), (2) Pareto-based 
method (two solutions), (3) multiple choice problem (two solutions), (4) HMMD (eight solutions). 



Table 21. Integration of composite solutions 



Method 


Resultant composite DAs 


Quality vector (HMMD) 


1 .Ideal-point 


= M 4 * L 2 * V 2 * U 3 * T 4 


(3; 2, 3,0) 


method 


S( = Mi ★ ii ★ V 2 * U 3 * Ti 


(3; 1,3,1) 




Si = Mi * L 4 * V 2 * f7 3 * T 4 


(3; 1,4,0) 




5| = M 4 * L 3 -kV 2 -kU 3 -kT 4 


(3; 1,4,0) 




Si = M h -kL 1 -kV 2 -kU 3 -kT 4 


(3; 1,4,0) 




Si = M 5 ★ L 2 * V 3 * U 2 * T 4 


(3; 1,2, 2) 




Si = M 5 * L 4 * K 2 * U 3 * T 4 


(3; 1,4,0) 


2. Pareto-based 


Sf = M 4 * L 2 * Vi * *7 5 * 


(3; 3, 1,1) 


MA 


S£ = M 4 *L 2 *V 2 *U 3 *T 4 


(3; 2, 3,0) 


3. Multiple choice 


Sf = M 4 *Li*V r 6 *Lf 3 *T 1 


(0;2,1,2) 


problem 


= M 4 ★ ii ★ y 6 ★ Z7i * Ti 


(0;3,0,2) 


4.HMMD 


S*f = M 4 *L2*V r B *£/'i*T 8 


(2; 5, 0,0) 




= M 4 ★ L 2 ★ V 5 ★ C/ 4 ★ T 2 


(2; 5, 0,0) 




S 1 ^ = M 4 * L 2 * Vi ★ C/ 5 ★ Ti 


(3; 3, 1,1) 




5f = M 4 ★ L 4 ★ y 5 ★ C/i * T 5 


(2; 5, 0,0) 




= M 4 * L 4 * y 5 * u 4 * T 2 


(3; 3, 1,0) 




= M 4 *L i *V 1 *U 5 * Ti 


(3; 3, 1,1) 




5f = M 4 *L 2 *V 2 *U 3 *T 4 


(3; 2, 3,0) 




Sf = M 4 * L 4 * V 2 * C/ 3 * T 4 


(3; 2, 3,0) 



Now let us discuss the obtained solutions (Table 21): 

1. In the case of the first three methods (MA, ideal point method, and Pareto-based method), compat- 
ibility estimates at level 3 were used to combine solutions. Thus cardinality of the combinatorial space 
of admissible solutions was decreased (for the examples) and the resultant solution set does not involve 
solutions with compatibility estimates at level 2 and 1. In the other case, cardinality of the admissible 
solution set can be very high. High cardinality of the admissible solution set will lead to very high 
computational complexity (MA, ideal point method, Pareto-based method) and participation of domain 
expert (s) at the first method stage (i.e., generation of admissible solutions) will not be possible. 

2. In the case of MA, a sufficiently large and rich set of admissible solutions was obtained: 45. Note, 
this solution set covers solutions sets for other methods (i.e., ideal-point method, Pareto-based method, 
HMMD). At the same time, the problem is: to analyze this large solution set. 

3. In the case of ideal point method, only solution Si belongs to the set of Pareto-efficient solutions. 
The set of considered solutions {S[ , S 2 , Si, 5|, Si, Si}, which are close to the above-mentioned solution, 
is not sufficiently good by elements. At the same time, some good solutions are lost, for example: S 3 , 

qH qH qH 

4- In the case of Pareto-based method, many good solutions are lost, for example: S$ , S§ , Sg , etc. 

5. In the case of multiple choice problem, compatibility estimates are not examined. As a result, all 
obtained solutions are inadmissible. It can be reasonable to extend this kind of optimization models by 
additional logical constraints which will formalize the compatibility requirements. But it may lead to 
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complicated models. 

6. In the case of HMMD, the set of solutions is sufficiently rich and not very large at the same time 
(eight solutions). 

Table 22 contains an additional qualitative author's comparison of the examined methods. Here com- 
putational complexity is depended on enumerative computing and analysis of all admissible combinatorial 
solutions (i.e., admissible combinations). In the case of HMMD, the usage of hierarchical system struc- 
ture decreases complexity of the computing process. In the case of Pareto-based MA, an analysis of 
Pareto-cfficient solutions will required additional enumerative computing. Finally, column "Usefulness 
for expert(s)" (Table 22) corresponds to the following: (i) possibility to include the domain(s) expert(s) 
or/and decision maker(s) into the solving process (i.e., to include cognitive man-machine procedures into 
the design framework), (ii) understandability of the used design method to domain(s) expert (s) and/or 
decision maker (s). 



Table 22. Qualitative comparison of used methods 



Method 


Computational 
complexity 


Taking into 

account 
compatibility 


Usefulness for selection 
of the best solutions 


Usefulness 
for 

expert (s) 


l.MA 


High 


Yes, binary 


Hard 


Hard 


2. Ideal-point method 


High 


Yes, binary 


Easy 


Good 


4.Pareto-based MA 


High 


Yes, binary 


Medium, analysis of 


Good 








Pareto-efficient solutions 




4. Multiple choice problem 


Low/Medium 


None 


Easy 


Medium 


5.HMMD 


Low/Medium 


Yes, ordinal 


Easy 


Good 



Generally, the selection of the certain kind of the morphological method for a designed system has to 
be based on the following: (a) a type of the examined system class (structure, complexity of component 
interaction, etc.); (b) structure and complexity of the examined representative of the system class; (c) 
existence of an experienced design team; (d) possibility to implement some assessment procedures (for 
assessment of DAs and/or compatibility); (e) possibility to use computational recourses (e.g., computing 
environment, power software, computing personnel), and (f) possibility to use qualified domain(s) experts 
and/or decision makers. 

5. Towards Other Approaches 

Generally, hierarchical design approaches are often based on a hierarchical model of the designed 
system and 'Bottom-Up' framework (Fig. 15). The list of some hierarchical design approaches, which 
are close to MA-based approaches and based on the framework above, is the following: (1) hierarchical 
design frameworks (e.g., [30], [69]); (2) structural synthesis of technical systems based on MA, cluster 
analysis, and parametric optimization [55]; (3) HTN (hierarchical task network) planning (e.g., [TS]); (4) 
hierarchical decision making in design and manufacturing (e.g., [5J, [6], [8], [23], [31]); and (5) linguistic 
geometry approach (e.g., [72]). 

Here it is reasonable to point out some non-linear programming models which are targeted to modular 
system design as well. First, modular design of series and series-parallel information processing from the 
viewpoint of reliable software design while taking into account a total budget (i.e., multi- version software 
design) was investigated in ([I], [2], [9]). The authors suggested several generalizations of knapsack 
problem with non-linear objective function. Thus, the following kind of the optimization model for 
reliable modular software design can be examined (a basic case) [9]: 

m qi m qi qi 

max ] [ (1 - JJ(1 - PijXij)) s.t. ^2 ^2 ,/ <-' •''<..' - & ' x v - 1 V * = * 

, 777, X{j €{0,1}, 

i=l j=l i=l j = l j = l 

where py is a reliability estimate of software module version (i.e., version j for module £), dij is 

a cost of software module version Fig. 16 illustrates the design problem above. Evidently, the 

obtained models are complicated ones and heuristics or enumerative techniques are used for the solving 
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process (pQ, [2], [9]). In [76], the problems above are considered regarding the usage of multi-objective 
genetic algorithms. Second, design problems in chemical engineering systems require often examination of 
integer and continuous variables at the same time and, as a result, non-linear mixed integer programming 
models are formulated and used (e.g., Q7], |21j). 



1. Composition/combination 
of the best DAs 

2. Analysis of resultant 
composite DAs 



T 



Composition/combination 
of the best DAs 



Composition/combination 
of the best DAs 



Selection of 
the best DAs 



I 



Selection of 
the best DAs 



Generation/ 
definition of DAs 



I 



Generation / 
definition of DAs 



Selection of 
the best DAs 



I 



Selection of 
the best DAs 



Generation/ 
definition of DAs 



I 



Generation/ 
definition of DAs 



Fig. 15. 'Bottom-Up' design scheme 




(a) series scheme (b) series-parallel scheme 

Fig. 16. Modular design of series or series/parallel system 



In addition, it is reasonable to point out constraint-based approaches (e.g., [TBJ, [55], [75]) including 
composite constraint satisfaction problems and Al-based solving approaches (e.g., [55], [7T]). 



6. Conclusions 

In the article, several MA-based system design approaches were described. Generally, it can be very 
useful and prospective to extend studies in the examination and usage of the MA based approaches in 
engineering, computer science, and management. For example, the following significant applied domains 
may be pointed out: (i) usage of morphological methods in allocation (layout, positioning) problems (e.g., 
[33) . |39j ) ; (ii) usage of morphological methods in combinatorial evolution and forecasting of modular 
systems (e.g., [57], [53]). The future research directions can include the following: 

1. continuation of the analysis, evaluation, comparison of MA-based system design methods; 

2. consideration of uncertainty in all modifications of MA; 

3. extension of "method of closeness to ideal point" while taking into account the following: (i) a set 
of ideal points, (ii) various kinds of proximity (e.g., functions, vector functions); 

4- analysis, investigation, and modification of morphological methods based on multiple choice problem 
and its generalizations including special constraints for system elements compatibility; 

5. design and investigation of special computer-aided systems based on morphological approaches; 

6. investigation of special versions of morphological approaches which involve experts into a solving 
process (i.e., interactive approaches); 

7. investigation of dynamical versions of morphological approaches while taking into account changes 
of system requirements; 
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8. usage of morphological system design methods for integration of heterogeneous networks; 

9. usage of morphological system design methods in embedded systems for configuration and reconfig- 
uration (including online mode) of hardware and/or software; 

10. generation of engineering benchmarks for evaluation and analysis of MA-based system design 
methods; and 

11. usage of MA and its modifications in engineering, IT/CS, and mathematical education (e.g., |37j . 
ED], 03], 09], [55, [52], [55], [56]). 
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